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Abstract 

Wavefunction at the origin with the incorporation of relativistic 
effect leads to singularity in a specific potential model. To regularise 
the wavefunction, we introduce a short distance scale here and use it 
to estimate masses and decay constants of and B s mesons within 
the QCD potential model. These values are then used to compute the 
oscillation frequency Arris of B& and B s mesons. The values are found 
to be in good agreement with experiment and other theoretical values. 
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The investigation of weak decays of mesons composed of a heavy quark 
and antiquark gives a very important insight in the heavy quark dynamics. 
The study of mixing and decay constants of the B meson, provides us useful 
information about the dynamics of quark and gluons at the hadronic scale. 
The weak eigenstates of neutral mesons are different to their mass eigenstates. 
This leads to the phenomenon of mixing whereby neutral mesons oscillate 
between their matter and antimatter state. This was observed first in the 
Kaon sector and subse-quently in Bd and B s mesons. The mass difference 
ArriB is a measure of the frequency of the change from a B into a B and so 
called the oscillation frequency. The decay constants of heavy mesons, one 
of the input parameters for oscillation frequency are crucial for interpreting 
data on particle-antiparticle mixing in the neutral B meson system, and for 
anticipating and interpreting new signatures for CP violation. 
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If the CKM element is well known from other measurements, then the pseu- 
doscalar decay constant f p can be well measured. If, on the otherhand, the 
CKM element is less well or poorly measured, having theoretical input on fp 
can allow a determination of the CKM element [1]. A measurement of the 
decay constant fp is difficult, since B + — > l + T*i is cabibo-suppressed in the 
Standard Model. Hence f Bq has to be provided from theory. 
Here in this paper we calculate pseudoscalar masses M Bq and pseudoscalar 
decay constants f Bq to compute the oscillation frequency ArriB q ,q — d,s 
within the frame work of a potential model [2,3,4]. To incorporate relativis- 
tic effect, the necessity of a short distance scale, in analogy to QED is also 
pointed out here. 

For the light heavy flavour bound system of qQ or qQ the hamiltonian 
can be written as 

H = ~ + V(r) (1) 
V (r) is the spin-independent quark-antiquark potential. 

V (r) = V coul (r) + V conf (r) (2) 

where V cou i represents coulombic part of the potential and V con f represents 
the confining potential. The vector and scalar confining potentials in the non 
relativistic limit reduce to [5,6] 

K n ,(r) = (l-e 1 )(br + C) (3) 

V: onf {r) = 6i (br + C) (4) 



reproducing 
and 



V conf (r) = V: onf (r) + V: onf (r) = br + C (5) 



V cou i(r) = -a c /r (6) 

where a c = ^a s , a s being the strong running coupling constant, e\ is the 
mixing co-efficient and b, C are the potential parameter as is used in our 
previous work [3, 4]. 
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Considering the linear part of the potential as perturbation, coulombic 
part as parent and then by using dalgarno method, the wavefunction in the 
ground state is obtained as[2,3,4 ] 



N'= . 2 " (8) 

^/(2*r (3 - 2c) C" 2 - I^r (5 - 2c) C" + ^/i 2 6 2 agr (7 - 2c)) 

C" = 1 + cAoV^g (9) 

(10) 



mi + rrij 
3' 



a = ( ^fjux^) (11) 



c = 1_ V 1_ v3 a v (12) 

Here is the undetermined factor appearing in the series solution of the 
Schrodinger equation, \x is the reduced mass and m,i,m,j are the costituent 

quarks masses .The term is due to the relativistic effects. The strong- 

running coupling constant appeared in the potential V(r) in turn is related 
to the quark mass parameter as [5] 

®s M = TV (13) 



(11-^) In 



where, n/ is the number of flavours, \i is renormalistion scale related to the 
constituent quark masses as a = 2 mimj and A is the QCD scale which is 
taken as 0.200 GeV here. Mb is the background mass related to the con- 
finement term of the potential as Mb = 2.24 x 6 1 / 2 = 0.95GeV '. In this 
calculation we have taken rif = 3 as in ref.[5,7]and the input mass param- 
eters as nid = 0.36GeV, m s = 0.46GeV and m& = 4.95GeV. With these 
values we calculate a s = 0.40 for B d mesons and a s = 0.37 for B s mesons. 
However we have also computed the coresponding results for rif =4 also and 
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listed in the tables. 

One of the problems faced in the model to study masses and decay con- 
stants is the incorporation of relativistic effect in the wave function at the 
origin, since a singularity develops at r = (eg. 7). However singularities at 
r = 0, in relativistic and nonrelativistic approach of quark model [5,6,8] is 
not new and different regularisation of singularities has been discussed. 
Here, in this work we use another way to regularise the wave function at the 
origin which have the quantum mechanical origin in QED. It is well known 
that relativistic wave function of the hydrogen atom too has such singulari- 
ties. However such effect is noticeable only for a tiny region[H], 

I 1 \ _ 2 _ 16300 

Imzar < e ^ 1 —i-><e ~ 10 ^ (14) 

Where z is the atomic number, m is the reduced mass of the hydrogen atom, a 
is the electromagnetic coupling constant and 7 = \J\ — z 2 a 2 . Using such hy- 
drogen like properties in QCD, m, a and 1 — 7 are to be replaced by /1, -^a s 



and e respectively. Here a s is strong coupling constants, e = 




and {mzar) 1 ' changes to y^A , leading to a cut off parameter r upto 
which the model can be extrapolated (r > tq). 

In analogy to the QED calculation (eg. 14), using the typical length scale for 
the relativistic term (jA < ^ , we get the cut off parameter 

r ~a e"^. (15) 

Unlike QED, it is flavour dependent depending on the flavours of the quark 
masses(eg.ll). 

The decay constants of mesons are important parameters in the study of 
leptonic or non-leptonic weak decay processes and in the neutrali?— B mixing 
process. In the nonrelativistic limit, the decay constant can be expressed 
through the ground state wavefunction at the origin ip p (0) by the Van-Royen- 
Weisskopf-formula [10]. Though most of the models predict the mesonic mass 
spectrum succesfuly, there are disagreements in the predictions of there decay 
constants. So we reexamine the predictions of the decay constants with a 
new short distance scale. 

We consider the nonrelativistic expression for f p as [10,11] 
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where M p is the pseudoscalar mass of mesons and is calculated by using 
the relation[12] 

M p = m t + m,-p^\^(0)\ 2 (17) 
6rriimj 

The values of ro and coresponding wavefunction at the origin for Bd and B s 



Table 1: Values of r and M p for B d and B s mesons 

i> (0) | in GeV z ' 2 values of r in Gev' 1 M B in GeV 



mesons 



B d 0.141 0.009 5.273 [our work with n f = 3] 

0.162 0.021 5.256[our work with n f = 4] 

5.279 [13] 
5.285 [14] 

5.279 [15] 

B s 0.178 0.002 5.370 [our work with n f = 3] 

0.207 0.007 5.349 [our work with n f = 4] 

5.369[13] 
5.373[14] 
5.375[15] 



are clculated to study pseudoscalar masses and are shown in table 1. Using 
eq.16, we compute the decay constants for Bd and B s mesons and put in 
table 2. The reults are found to be well with the experimental values. 

The neutral Bd and B s mesons can mix with their antiparticles by means 
of box diagram involving exchange of a pair of W bosons and intermedi- 
ate u, c, t quarks leading to oscillations between the mass eigenstates. This 
mass oscillation is parametrized as the oscillation frequency or mixing mass 
parameter(Am) and is given by [16,17] 

Am B = \ qJB " g (x t ) Vt\V t * q V tb \ 2 B, q = d, s (18) 
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where r] t is the gluonic correction to the oscillation and is taken 0.55 as in 
ref.[16], the last factor B is the bag parameter which represents the correction 
to the vacuum insertion and is taken as 1.34 [16]. The function g (x t ) is given 
by[18] 

, < 1 9 3 3^Cj. , . 

9{x ' ) = l + W^)-w^f~W^ (19) 

Here, x t = The values m t (17AGeV), m w (80.403Ge\/),and the CKM 

W 

matrix elements \V tb \ (1), |14i| (7.4 x 10~ 3 ),|Vt s | (40.6 x 10~ 3 ) are taken from 
the particle data group [13]. 

We use the estimsted values of masses and decay constants to compute the 
mixing mass parameters and put in the table 2. The decay constants are 
found to be comparable with the other theoretical values and mixing mass 
parameters are found to be in good agreement with the experimental values. 



Table 2: Decay constant and oscillation freq. of B mesons 



mesons /fi(i n GeV) Am#(in ps r ) 



Bd 0.213[our work with rif = 


3] 


0.55[our work with rif ■ 


= 3] 


0.246 [our work with rif — 


4] 


0.74 [our work with rif ■ 


= 4] 


0.189 [22] 




0.50[13] 




0.196(29)[23] 




0.547[20] 




0.190(7)±f [22] 




0.515[21] 




0.210(9) [25] 








0.216(9)(19)(6)[26,27] 








B s 0.265 [ our work with rif — 


= 3] 


17. 34 [our work with rif 


= 3] 


0.311 [our work with rif = 


4] 


23.88[our work with rif 


= 4] 


0.218 [22] 




17[13 ] 




0.216[23] 




17.77(19] 




0.217(6)lg[24] 









0.244(21) [25] 
0.259(32) [26,27] 



To conclude, in the present work, we have incorporated the relativistic 
effects to the wave function at the origin in a QCD potential model. We have 
introduced a short distance scale r in analogy to QED. Unlike QED, in QCD 
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such short distance scale become flavour dependent. As expected, its mag- 
nitude is far larger than its QED counterpart but still far smaller than the 
measure of finite size of hadrons or its constituents. It is however well within 
the reach of LHC[26] where distance down to a scale as short as 5 x 10~ 20 m 
-lCT 21 m will be explored. Theoretically this short distance scale r can be 
roughly associated with the ultraviolate regularisation scale of QCD. 

We have found from the calculation that generally the masses are not 
so sensitive to running coupling constants since the heavy quark constituent 
masses dominent here, but the decay constants are very much sensitive to 
the running coupling constants and therefore a slight change in a s as well as 
r deviates the results significantly. 

In table 2 we confront our results for pseodoscalar decay constants and os- 
cillation frequency(Am#) of the Bd and B s mesons with the recent pre- 
dictions of quark model [22], three flavour lattice QCD [27,28], Collabora- 
tion^, 20, 21] and available experimental data[13]. In our calculation we have 
found ^=1.24 which is in accordance with the lattice result ^-=1.20(3) (1) 
[26,27]. The results in table 1 and table 2 for Bd and B s mesons with the 
short distance scale is found to be in good agreement between all presented 
theoretical predictions as well as experimental data for rif = 3 and the cor- 
responding results with rif — 4 are found to overshoot the theoretical pre- 
dictions and experimental data. 
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